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In this contribution, we study situations in which nanoparticles in 
a fluid are strongly heated, generating high heat fluxes. This sit- 
uation is relevant to experiments in which a fluid is locally heated 
using selective absorption of radiation by solid particles. We first 
study this situation for different types of molecular interactions, us- 
ing models for gold particles suspended in octane and in water. As 
already reported in experiments, very high heat fluxes and temper- 
ature elevations (leading eventually to particle destruction) can be 
observed in such situations. We show that a very simple modeling 
based on Lennard-Jones interactions captures the essential features 
of such experiments, and that the results for various liquids can be 
mapped onto the Lennard-Jones case, provided a physically justified 
(corresponding state) choice of parameters is made. Physically, the 
possibility of sustaining very high heat fluxes is related to the strong 
curvature of the interface that inhibits the formation of an insulating 
vapor film. 
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Introduction 

Sub-micron scale heat transfer is attracting a growing inter- 
est, motivated by both fundamental and technological points 
of view. In fluids, considerable attention has been devoted to 
the so called nanofluids jT] [2] , in which nanoparticles in dilute 
suspension appear to modify both bulk heat transfer and crit- 
ical heat fluxes. While the former effect can presumably be 
understood in terms of particle aggregation [3 , 4 , the latter 
is still poorly understood. 

More generally, heat transfer from nanoparticles or nanos- 
tructures to a fluid environment is a subject of active research, 
stimulated by the development of experimental techniques 
such as time resolved optical absorption or reflectivity, or pho- 
tothermal correlation spectroscopy |5j. Applications include 
e.g. the enhancement of cooling from structured surfaces, lo- 
cal heating of fluids by selective absorption from nanoparti- 
cles, with possible biomedical hyperthermia uses [6l [7]. Re- 
cent experiments demonstrated the possibility of reaching very 
high local temperatures using laser heating of nanoparticles 
[3 El Unj? even reaching the melting point of gold particles 
suspended in water. From a conceptual point of view, such ex- 
periments raise many interesting questions compared to usual, 
macroscopic heat transfer experiments. How are the phase di- 
agram and heat transfer equations modified at small scales ? 
How relevant is the presence of interfacial resistances and how 
do they change with temperature ? 

The case of nanofluids [11 is a good illustration of the role 
that can be played by molecular simulation in the interpreta- 
tion of such complex situations. While many interpretations 
have been proposed to explain the reported experimental re- 
sults, it is only simulation of simple models that has been able 
to disprove some of these interpretations and to demonstrate 
the validity of the alternative, aggregation scenario. Inter- 
estingly, the use of complex models with accurate interaction 
force fields is not, in general, needed to answer the basic quali- 
tative questions raised by such novel experimental approaches. 



In this manuscript, we use molecular simulation to study 
the heat transfer from solid nanoparticles to a surrounding 
fluid under extreme conditions (high heat fluxes) using both 
realistic and simplified molecular models. We show that, in 
agreement with experiments, the temperature of nanometer 
sized particles can be elevated considerably without inducing 
bubble nucleation in the fluid. This feature is contrasted with 
the situation for flat surfaces, at which an instability leading to 
the formation of an insulating vapor layer takes place at much 
smaller heat fluxes and temperatures [12| [T3] . Using a com- 
parison between a "realistic" description of gold-octane and 
gold-water systems and of a simplified Lennard-Jones model, 
we show, based on a "corresponding state" analysis, that the 
features observed are quite universal. A simple mapping using 
critical temperature, interfacial and heat conductivity prop- 
erties of the particle/solvent pair allows one to reproduce ac- 
curately the behavior of different systems. 



Gold nanoparticles in octane 

To establish a connection with experimental reality, we began 
by using molecular dynamics to simulate heat transfer from 
a model gold nanoparticle into octane solvent. The selection 
of a relatively simple organic solvent in this first part, rather 
than water, is motivated by the fact that molecular models are 
quite accurate in predicting liquid-vapor phase diagrams for 
alkanes. Our model consists of a nanoparticle, with an average 
radius of about 1.3 nm, containing 494 gold atoms arranged 
on a FCC lattice with density 19.5 g/cm^. The nanoparticle 
is immersed in liquid octane containing 2721 octane molecules 
(21768 united carbon- hydrogen atoms) and placed in a cubic 
simulation cell with periodic boundary conditions. The inter- 
actions between united atoms forming octane molecules are 
described by the Amber force-field |M with all non-bonded 
interaction energy calculated according to 6-12 Lennard-Jones 
(LJ). The interaction between gold atoms was also described 
by a 6-12 LJ potential, fitted to reproduce bulk Au density 
(19.5 g/cm^) and the melting point (1310 K). Finally, the 
interaction potential between octane united atoms and gold 
atoms was adapted from Ref. [15 . All simulations were car- 
ried out at constant pressure of 1 atmosphere and with the 
integration time step of 2 fs. In the equilibration stage of 
the simulation, a global thermostat is used to maintain the 
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overall temperature at 300 K. The equilibration stage takes 
approximately 1 ns, and under 1 atm, the equilibrated system 
is contained in a cubic simulation box with the edge length of 
about 90 A. The nanoparticle is initially placed in the center 
of the simulation cell, but is allowed to freely diffuse during 
the simulations. . The equilibrium density of the model oc- 
tane fluid is equal to 0.71 g cm~^, which compares well with 
the experimental density of 0.7025 g cm~^. To study the heat 
flow from the nanoparticle to the solvent, the nanoparticle 
was heated, with a constant heating power in the range of 
100-1000 nW, by rescaling of the atomic velocities every time 
step. The liquid octane in the periphery of the system, at dis- 
tances more than 40 A from the nanoparticle center (taking 
into account a possible diffusive motion), was maintained at 
300 K, thus providing the heat sink. Up to heating powers of 
700 nW, after a transient of about 100 ps, a steady state is es- 
tablished. In the steady state, we collect time averages (over 
Ins) of density and temperature profiles obtained for spher- 
ical shells concentric with the nanoparticle center and with 
a thickness of 2A. For heating powers larger than 700 nW, 
the system was unstable and its behavior will be described 
below. 
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Fig. 1. (Color online) Steady state temperature profiles for octane-gold model 
system at three heating power levels (symbols) and fits with the continuum theory 
prediction T(r) = A -\- B/r (solid lines) 

Steady state temperature profiles for P = 200, 400 and 
600 nW heating powers are presented in Fig. [T] The tem- 
perature profiles have several noteworthy features. First, the 
temperature of the nanoparticle is more or less uniform; this 
is an effect of the relatively high thermal conductivity of crys- 
talline solid, as compared to liquid. We mention here, that 
the electronic contribution to the conductivity, which provides 
the dominant mechanism of heat conduction within metallic 
nanoparticles, is not accounted for in our description. The 
essential point, however, is that the conductivity of the solid 
is much higher than that of the liquid, independently of the 
precise mechanism involved. Second, the temperature in the 
liquid follows the solution of the continuum heat flow prob- 
lem with spherical symmetry (see solid lines in Fig. [1]). In 
the steady state the temperature profile is described by the 
solution of the Laplace equation in the form: 

T(r) = A + B/r [1] 

Near the particle-liquid interface, the temperature profile de- 
viates from the formula given by the Eq. ([T]), particularly for 
larger heating powers involved. This deviation is likely due 
to the non-uniform thermal transport properties of the liquid, 
since Eq. [T] is valid under the assumption that the thermal 



conductivity is constant. Very importantly, there is a large 
temperature drop, AT, at the nanoparticle- liquid interface 
which is a manifestation of the interfacial thermal resistance. 
Such resistance is caused by the mismatch of thermal prop- 
erties between the solid and liquid components, and is also 
affected by the strength of the interfacial bonding. The in- 
terfacial thermal conductance, C, can be quantified via the 
relationship: 



jQ = GAT [2] 



Where jq is the heat flux across the interface and AT is the 
discontinuous temperature at the interface (see temperature 
profiles in Fig. [T]). Fig. [2] (top panel) shows the relation- 
ship between the heating power P and the temperature drop 
at the octane-liquid interface. At lower heating powers the 
heat flux is proportional to the heating power, indicating con- 
stant value of the interfacial conductance. However, above the 
300nW heating power, the increase in the temperature drop 
becomes steeper, indicating increasing interfacial thermal re- 
sistance (see Fig. [2]). 
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Fig. 2. (Color online) Temperature drop at the interface for octane-gold model 
system (top panel points) and calculated interface conductance (bottom panel points) 
as a function of heating power. Lines represent linear response regime. 

The calculated interfacial conductance as a function of 
P is shown in Fig. [2] (bottom panel). At small heating 
powers (temperature drops), the value of the interfacial con- 
ductance is about 100 MW/m^/K. This value is similar to 
those obtained in an experiment on gold nanoparticle-water 
dispersions [16]. With increasing heat power (temperature 
drop), the interfacial thermal conductance decreases from 100 
to about 50 MW/m^/K at P = 700nW. To gain an insight 
into the structural origin of the behavior, we show in Fig. [3] 
the octane density profiles corresponding to temperature pro- 
files from Fig. [T] As the temperature of the nanoparticle and 
the adjacent liquid increases, there is a visible decrease of liq- 
uid density adjacent to the solid surface (see Fig. [3)). This 
increase of the molecular distance between liquid molecules 
and solid atoms is likely responsible for the decrease of the 
interfacial thermal conductance. 
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scope of this article, and we leave a complete study for future 
work. 
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Fig. 3. (Color online) Steady state octane density profiles for octane-gold model 
system at three heating power levels. 

Gold nanoparticles in water 

In order to explore a different type of bonding for the liquid, 
and also to allow a connection to recent experiments [5 , we 
explore in this section a system made of a gold nanoparticle 
similar to the one described in the previous section, solvated in 
water modeled using the standard SPC/E model [17]. In order 
to allow comparisons with the Lennard- Jones calculations de- 
scribed below, we studied water pressurized at a pressure of 80 
bars, and at a temperature of 450 K. These parameters were 
chosen so that the distance to the critical point is large, but 
still within a range that allows comparisons between different 
models. The system was made of 10000 water molecules. 

The main issue here is the choice of the model for the gold- 
water interaction. Experimentally, the results reported for the 
wetting of water on gold are quite scattered, with in general a 
large contact angle hysteresis. Results obtained under UHV 
conditions [18 , 19 report a low contact angle (smaller than 30 
degrees, or close to a wetting situation). On the other hand, 
the force fields that exist in the literature, and have been 
based on density functional theory DFT calculations p!2l |2Q] . 
yield higher contact angles 

In this study we make the choice of strengthening the at- 
tractive terms in this effective potentials, in order to obtain 
a contact angle (estimated from a simple calculation of the 
Lennard- Jones contributions to the surface tensions [21 , and 
of a the actual value of the SPC/E surface tension at 300K) 
of the order of 25 degrees, consistent with experiments. As a 
result, the gold- water interaction is written in the form of a 
standard 6-12 potential like in 120] . with the following param- 
eters: toiAu — 0.59 kcal.moP and gqiau — 0.36 nm, while 
the hydrogen atoms do not interact with the gold atoms. The 
system NP + water is first equilibrated during 100 000 time 
steps which represent a physical time of 200 ps. Then the 
nanoparticle is heated up at a constant power while the water 
molecules at a distance 20 Afrom the nanoparticle are ther- 
mostated at 450K. In all the following, we will restrict our- 
selves to moderate heating powers (smaller than 700 nW), 
since for larger heating intensities, we have observed non- 
stationary effects in the heat transfer process. While these 
effets are interesting in themselves, their study is out of the 
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Fig. 4. (Color online) Temperature profiles across the water-gold nanoparticle 
interface at three heating strengths. 

Temperature profiles of gold NP immersed in water, 
heated at various powers, are shown in figure 2] Clearly, tem- 
perature profiles are flatter in water than in octane, if the 
comparison is made at the same value of the heating power. 
This is due to the about 5 times larger conductivity of water 
compared to octane. On the other hand, interfacial temper- 
ature jumps AT are smaller in the water/gold case. Note 
that it is essential to do the comparison at a given value of 
the heating power, not at a given value of the nanoparticle 
temperature. For instance, for the nanoparticle heated up at 
400 nW in octane, we have measured AT = 220 K, while it 
is a factor of 2 less if the nanoparticle is immersed in water. 
Consequently, the water/gold interface has a larger interfa- 
cial conductance than the octane/gold system. The value we 
have measured varies from G = 170 MW/m^K to G = 150 
MW/mVK over the range of heating power investigated, a 
variation that is smaller than in the case of octane. This 
trend is consistent with the recent finding that the interfacial 
conductance increases with work of adhesion, which is higher 
for the gold- water case than for the gold-octane case. 22 The 
Kapitsa length Ik — ^/G^ where A is the thermal conductiv- 
ity, is of the order of 3.4nm in this case and of Inm in the 
gold/octane case. 



Melting of the nanoparticles 

At high enough temperatures, experiments have illustrated 
the possibility of particle melting [9] within the fluid. We also 
explored briefly this issue in our simulations of gold particle in 
octane. By inspection, we observed that up to P = 500 nW, 
the nanoparticle structure remained crystalline. However, at 
P = 600 nW, the crystalline order of the nanoparticle is lost. 
At P = 700 nW, we observed that atoms from the nanopar- 
ticle surface are gradually "evaporated" into the solution. At 
a later stage they recombine into small Au clusters (Fig. [5]). 
We note that all these processes occur without formation of 
a liquid vapor interface, which, as discussed below, is caused 
by very large Laplace pressure. In fact, using the surface ten- 
sion of octane at room temperature 7 = 21.8 10~^ N/m and a 
bubble radius of i?o=2 nm, one obtains the Laplace pressure 
of Pl=2j / Ro ^ 200 atm. This value is much larger than the 
critical octane pressure of 25.5 atm. These results illustrate 
an exciting possibility of decomposition of metal nanoparticles 
into metal atoms or small clusters, without explosive evapora- 
tion or thermal damage on the embedding medium. Interest- 
ingly, we have not observed such fragmentation for nanoparti- 
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cles immersed in water. Although the crystalline order is lost, 
the gold nanoparticle keeps its integrity. This is probably due 
to the higher interfacial energy of the gold/water interface 
compared to gold/octane. 
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stat). The temperature To is below the boiling temperature, 
that we found to be Tb ^ 0.8, using independent simulations of 
a liquid/vapor interface, under the pressure Po we are working 
at. After 100000 time steps of equilibration, the nanoparticle 
is heated up at different temperatures Tp > Th by rescaling 
the velocities of the solid particles at each time step, while 
the whole system is kept at the constant pressure Po using 
a NPH barostat. The fluid beyond a distance lOcr from the 
particle surface is thermostatted at To = 0.75, again using ve- 
locity rescaling. Temperature, density and pressure fields have 
been obtained by averaging the corresponding quantities dur- 
ing 10000 time steps in nanoparticle centered spherical shells 
of width ~ 0.15(7, after a steady state is reached. Finally, 
we calculate the heat flux density flowing through the solid 
particle, by measuring the power supply needed to keep the 
nanoparticle at the target temperature Tp. 
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Fig. 5. (Color online) Snapshots of octane-gold model system at heating powers 
of nW (a), 500 nW (b), 700 nW (c) and 1000 nW (d). 

Lennard-Jones model 

To emphasize the generality of the scenario described above, 
we briefly recall here the results obtained for a generic model 
of a simple, Lennard Jones monoatomic fluid in contact with 
a heated solid particle p^. The system is made of a FCC 
solid particle formed of 555 atoms, immersed in a fluid of 
23000 atoms, all the atoms interacting through a Lennard- 
Jones potential Va/3(r) = 4e((cr/r)"'^^ — Ccxpi^cj /r)^) where a^[3 
refers to solid or liquid atoms. The potential has a cut-off 
radius 2.5cr where a is the diameter of the atoms. The pa- 
rameters e and a are taken to be the same for both phases. 
The parameter Cap — 1 Ma — f3; Cap — cfs otherwise con- 
trols the wetting interaction between the fluid and the solid 
nanoparticle. In addition to the Lennard-Jones interactions, 
atoms inside the particles are connected to their neighbors 
with FENE springs [23^ V{r) = -0.5A;i?g In (l - (r/i?o)^) 
with k — 30e/cr^ and Rq — 1.5cr. This nearest neighbor 
bonding allows one to heat up the nanoparticle to rather high 
temperatures without observing the melting or fragmentation 
phenomena mentioned above.. In the following, we will con- 
centrate on results obtained for cfs — 1, which will be shown 
in section to be an appropriate value for a mapping between 
the LJ and gold-octane systems. The dependence on cfs is 
discussed in ref. [13 . All results in this section are given 
in standard Lennard-Jones units, e//cs, a and r = yjma'^ /e 
for temperature, length and time, respectively. Here m is the 
mass of the fluid atoms. 

We integrate the equations of motion using a velocity Ver- 
let algorithm with a time step dt — 0.005t. All the systems 
considered have been first equilibrated at a constant tempera- 
ture To = 0.75 under the constant pressure Po = 0.015 (using 
a Nose/Hoover temperature thermostat and pressure baro- 
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Fig. 6. (Color online) Steady state temperature field across the liquid/nanoparticle 
interface, obtained with the Lennard-Jones model (adapted from the results of ref. 
fl3l . The temperature has been normalised by the critical temperature Tc of the 
Lennard-Jones fluid. The position of the nanoparticle surface is indicated by dashed 
lines, r measuring the distance to the center of the nanoparticle. Solid curves cor- 
respond to fits by the continuum theory result T(r) = A/r + B. Note that for 
Tp = 1 and Tp = 1.5, the solid curves are almost indistinguishable from the sim- 
ulation data. For the highest temperatures, the temperature field inside the particle 
is not shown, in order to limit the amplitude of the scale of the vertical axis. 

Figure [6] displays steady state temperature profiles close 
to the nanoparticle surface, for different temperatures Tp of 
the nanoparticle. For low Tp, the temperature field in the 
liquid is practically indistinguishable from the form A/r -\- B 
(eq. [T]). The general behavior is strikingly similar to the one 
obtained for the gold-octane system. Inside the solid, the 
temperature is not uniform but slightly curved downwards, 
due to the finite conductivity of the nanoparticle. As in the 
gold/octane simulation, the temperature is discontinuous at 
the interface. For a given value of the nanoparticle tempera- 
ture Tp, the relative temperature jump AT is quite compara- 
ble to the gold-octane results, thus suggesting that the corre- 
sponding interfacial conductance is depends on the details of 
the nanoparticle/fluid interaction mostly trough macroscopic 
parameters such as the wettability ^241 1251 i26 .. Similar indi- 
cations can be found in Ref. [22 , where it was shown that 
the interfacial conductance of surfactant-water interfaces is 
directly proportional to the work of adhesion. For instance, 
for the nanoparticle heated up at Tp = 1.5, we have measured 
a temperature jump AT/Tc 0.4 which is quite compara- 
ble to the value reported in fig. [2] for the gold nanoparticle 
for the 400 nW heating power. Similarly, the value of the 
conductance obtained G = jq/AT ~ 0.6 is consistent with 
the gold/octane results displayed in fig. [21 if we assume that 
a value of C = 1 in our LJ units corresponds to an interfa- 
cial conductance on the order of 100 MW/K/m^. All these 
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results suggest that the interfacial conductance is a quantity 
which does not strongly depend on the details of the nanopar- 
ticle/fluid interaction but only on generic properties of an in- 
terface as e.g the wettability. This correspondance between 
the simplified Lennard- Jones system and the more realistic 
description of gold in octane will be explored further in the 
next section. 

We also mention briefly that very high flux situations can 
be explored here thanks to the "covalent" bonding introduced 
between the atoms of the solid particle. Upon increasing the 
temperature of the nanoparticle, deviations from the 1/r be- 
havior are clearly seen in fig. [G] when the local temperature 
exceeds the critical temperature Tc. Interestingly, the temper- 
ature profile steepens close to the nanoparticle surface, cor- 
responding to a decrease of the local effective conductivity.. 
The density profile [13] differs somewhat from its octane or 
water counterparts, with density oscillations in the case of 
the monoatomic LJ fluid. These oscillations are smeared out 
for a molecular liquid. The most important point is that, even 
far above the critical point, we do not observe a steep decrease 
of the Uquid density in the vicinity of the particle, but rather 
the appearance of a dilute liquid layer, with a density decreas- 
ing when the nanoparticle gets hotter. Note however that the 
density within this layer is still one order of magnitude larger 
than the vapor density at coexistence. Thus, even at temper- 
atures several times Tc, boiling of the surrounding fluid is not 
observed. 



Corresponding states analysis 

Obviously, the qualitative similarities between the observa- 
tions in the previous sections points to a rather generic sce- 
nario. In order to allow a more quantitative comparison, a 
mapping between the different systems is necessary. The map- 
ping we investigate will be based on the physical properties 
that dominate the problem under consideration, namely in- 
terfacial effects, liquid vapor coexistence and heat transfer in 
the Uquid phase. 

We start by a discussion of the gold-octane case, which 
in view of the rather similar interactions can be expected to 
be easily mapped onto a Lennard- Jones system. The first 
step is to tune the interaction coefficient between fluid and 
solid, cfs, to a value that is best suited to reproduce the 
properties of the gold octane interface. Generally speaking, 
the wetting properties of an interface can be related to the 
interaction potential Uij{r) through the interfacial work [21] 
^ij — — j^^^ ruij(r)dY where ro is a minimal radius of 
approach between two molecules, pi and pj denoting the num- 
ber densities of the interacting media. The wetting properties 
(equilibrium contact angle or spreading parameter) will be 
determined by the ratio r = Hps / Hff- In the gold octane 



system, this ratio is 



In the latter 



PCHI ^CH2/CH2 ^CH2/CH2 

expression, pAu =0.1 mol.cm"^ and pcH2 — 0-05 mol.cm"*^ 
are the gold and fluid number densities, eci^2/An =0-429 
kcal.mol"^; eci^2/Ci^2— ^.143 kcal.mol"^ are the octane united 
atoms/gold and octane/fluid interaction energies [15] , while 
c^Cif2/An=0.328 nm; crcH2/CH2—^-^^'^^ the radii of 

the corresponding interaction. Hence we find r ^ 2. For the 
monoatomic LJ model, r = cpsps/pF where the nanoparti- 
cle and liquid densities are, respectively, ps — 1.46cr~^ and 
pL — 0.76(7 To match the value of r, a value of parameter 
cfs — 1, as used in section , is therefore appropriate. . 

Now, we discuss the value of the unit of thermal flux used 
in the generic LJ model. To this end, we shall first determine 
the units of length, time and energy cr, r, and T correspond- 



ing to the generic model. Throughout, we denote by stars 
quantities expressed in LJ units where the units of length a, 
time T and energy e are all set equal to 1. Alternatively, 
we can determine the values of a, r and e by matching the 
thermal properties of the generic LJ model to the ones of 
the gold/octane model. By matching the critical temperature 
of the LJ model, T* = 1.08 to the octane critical tempera- 
ture Tc — 569 K, we obtain the unit of energy, e/ks — 527 
K. The values of the units of length and time are obtained 
by matching the values of the thermal conductivity A and 
thermal diffusivity Dth of octane cr^ = i^thkB and 

^ ^th'^B 

r = ^ ^ . The thermal conductivity of the monoatomic 

LJ fluid A"^ = 0.36 was measured using stationary heat trans- 
fer simulations, while we have used the value D^^ = 1 for 
the thermal diffusivity reported in [27]. Using the values 
of the thermal conductivity A = 0.1 W.m~^.K~^ and the 
thermal diffusivity Dth = 6.4 10"^ m^.s"^ at 400 K [28], 
we obtain a — 0.32 nm and r = 1.6 ps. The power of 1 
in Lennard Jones units corresponds then to = 5 nW, 

while a boundary conductance = 1 is equivalent to a real 
G = /cs/rcr^ = 88MW/K/m^ 
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Fig. 7. (Color online) Reduced temperature profiles across the gold/octane in- 
terface (filled symbols) compared to the profiles obtained with the LJ model (open 
symbols). Tc is the critical temperature of octane and of the LJ fluid respectively. 
For the gold/octane system, a = 0.32 nm as discussed in the text. Here the 
gold/octane system was thermostatted at T = 380 K. 



To illustrate the relevance of the mapping discussed, we 
have compared in fig. [7] the temperature profiles correspond- 
ing to the gold/octane interface and to the LJ model, in terms 
of the reduced temperature T /Tc. The distances have been 
rescaled here by the value of a = 0.32 nm discussed before. 
For the sake of the comparison, the gold/octane systems have 
been thermostatted at a higher temperature T — 380 K than 
before, since at the reduced temperature T = 300/569, the 
LJ fluid may cristallise. The agreement between the atomi- 
cally realistic model and the monoatomic LJ model is fairly 
good. For the lower heating strengths considered, the LJ fluid 
develops almost the same temperature profile away from the 
nanoparticle. In particular, the slopes of the temperature pro- 
file (flux) at the solid interfaces are quite comparable. Note 
however, that the temperature jump AT at the interface is 
larger in the case of the LJ model, a fact which can be at- 
tributed to the slightly smaller values of the interfacial con- 
ductances of the monoatomic model compared to the more 
realistic system. For higher heating powers, the LJ tempera- 
ture lies slightly above the octane curve. However, the fluxes 
at the solid interfaces are again quite comparable and the in- 
terfacial temperature jumps AT compare well, probably due 
to the rapid decrease of the gold/octane conductance with the 
supplied power. 
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We have carried out a similar analysis for the gold/water 
system. The values retained are A = 0.58 W/m/K and 
Dth = 1.3 10~^ m/s for the conductivity and thermal diffusiv- 
ity at T = 450 K, and Tc = 650 K for the critical temperature. 
This yields a = 0.22 nm, r = 0.384 ps and e/kB = 602 K. for 
the unit of length, time and energy, respectively. The values of 
the length and time units are smaller than their gold/octane 
counterpart, due to the larger thermal conductivity of water. 
As a consequence, the units of power P = 21.6 nW and in- 
terfacial conductance G = 537 MW/m^/K are larger than 
in the previous case. Because of the complexity of the in- 
teratomic potentials involved in the water/gold system, we 
have used a different approach than in the gold/octane case 
to adjust the soUd- liquid interaction parameter for the "corre- 
sponding" Lennard- Jones system. The reduced temperatures 
T/Tc of the nanoparticles where chosen to be identical, and 
the parameter c was chosen such that the power input is sim- 
ilar in both cases. This essentially amounts to matching the 
interfacial conductances of the two systems. The resulting 
value c = 0.78 leads to an interfacial work ratio r = 1.56. 
This is slightly different from the value r = 1.9 based on 
the gold/water interfacial work TvpAupo^o/Au^^o/Au/^ — 0-13 
kcal/mol/A^ and the water surface tension 7 ~ 40 mJ/m^ at 
450 K[17]. This difference is probably due to the more com- 
plex structure of the gold- water interface, which implies that 
a simple matching based on identifying the r parameters is 
not appropriate. 

When the matching is performed on the interfacial con- 
ductance, the agreement between the two temperature pro- 
files, shown in figlH is very good. . That this is the case may 
seem obvious, as the main macroscopic parameters have been 
matched. However, the fact that the agreement is obtained 
down to subnanometer scales, in a situation where heat fluxes 
are extremely strong, and can be transferred to different val- 
ues of the power input, is far from trivial. 
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Fig. 8. (Color online) Reduced temperature profiles across the gold/water interface 
(filled symbols) compared to the profiles obtained with the LJ model (open symbols). 
Here Tc = 650 K and a = 0.22 nm for water. 



In conclusion, the simple monoatomic LJ model, with an 
appropriate choice of the parameters, can be used le to re- 
produce quantitatively the features of heat transfer around a 
nanoparticle obtained for systems with more complex interac- 
tions. It constitutes a simple and efficient tool to explore heat 
transfer at the nanoscale, in th espirit of a coarse graining 
approach of such mesoscale phenomena. 



Conclusions 

We have explored the phenomenon of heat transfer in the 
vicinity of strongly heated nanoparticles, using molecular dy- 
namics simulations of atomically realistic models or of more 
coarse grained Lennard- Jones monoatomic fluids. The com- 
parison between the two approaches shows that, provided the 
mapping is carried out using the physically relevant proper- 
ties, they are quantitatively equivalent. The simulations re- 
veal that the fluid in the vicinity of the nanoparticles can sus- 
tain very high heat fluxes and large temperature differences 
without undergoing the type of drying instability that is ob- 
served on flat surfaces, and that temperatures much above the 
critical temperatures can be reached without observing phase 
coexistence . In the case of gold in octane, high heat fluxes 
and temperatures can result in a partial desintegration of the 
nanoparticle, while in the gold in water case they only result 
in melting of the particle for comparable heating powers.. 

The phenomena that are involved in such experiments are 
quite complex, with a combination of phase transition, inter- 
facial phenomena and transport phenomena, all taking place 
on the nanometer scale. It is therefore not trivial, that a 
simple coarse graining based on a matching of interface and 
thermal properties results in a correct description of the phe- 
nomenon. This implies that coarse grained methods should 
be appropriate for describing similar phenomena in nanostruc- 
tures of larger dimensions (e.g. aggregates of nanoparticles, or 
nanostructured surfaces) . The level of coarse graining and the 
gain in efficiency provided by the monoatomic LJ fluid could 
be even improved by extending to appropriately adapted ver- 
sions of other coarse graining approaches such as free energy 
models [29^ or appropriate versions of dissipative particle dy- 
namics 130 . 
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